Spin-Polarized Semiconductor Induced by Magnetic Impurities in Graphene 



O 
(N 

o 

O 

r— I 



I 

O 

o 



> 

rn 

o 
o 



X 



Maria Daghofer,^'^^LJ Nan Zheng, ^'■^ and Adriana Moreo^ 

'Department of Physics and Astronomy, University of Tennessee, Knoxville, TN 37996-1200 and 

Materials Science and Technology Division, ORNL, Oak Ridge, TN 37831-6032, USA 

^ IF W Dresden, P.O. Box 27 01 16, D-01171 Dresden, Germany 

^Department of Computer Science, The College of William and Mary, Williamsburg, VA 23187, USA 

(Dated: October 5, 2010) 

The effective magnetic coupling between magnetic impurities adsorbed on graphene, which is 
mediated by the itinerant graphene electrons, and its impact on the electrons' spectral density 
are studied. The magnetic interaction breaks the symmetry between the sublattices, leading to 
antiferromagnetic order, and a gap for the itinerant electrons develops. Random doping produces 
a semiconductor, but if all or most of the impurities are localized in the same sublattice the spin 
degeneracy can be lifted and a spin-polarized semiconductor induced. 

PACS numbers: 71.10.Fd, 72.25.Dc, 73.22.Pr 



I. INTRODUCTION AND MODEL 

Graphene has been the focus of active research in the 
last few years due to its peculiar transport properties'^— 
and its potential technological applications. For spin- 
tronic devices jS where the high electron mobility and tun- 
able electron filling make graphene an interesting candi- 
date, it is important to investigate magnetism and the 
possibility of spin-polarized conduction electrons. More- 
over, it would be desirable to transform graphene from a 
semi-metal into a semiconductor. In the present paper, 
we are going to show that localized magnetic impurities 
(Mis) provide a way to open a gap. 

Density Functional Theory (DFT) calculations sug- 
gest that decorating graphene-nanoribbons with metal 
ions can lead to opposite spin polarization in the valence 
and conduction bands. ^ This property is known as spin- 
polarized half-semiconductivity and it is very desirable for 
spintronics devices.— DFT also predicts that interfaces 
between Ni or Co layers and graphene or graphite could 
induce spin polarization^ because only the majority spin 
bands of the Co or Ni have an overlap with the graphene 
bands. We will show that Mis, which do not directly do- 
nate spin-polarized electrons to the graphene layer, can 
nevertheless polarize its electrons in the direction oppo- 
site to their own spin. 
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FIG. 1: (Color online) Cartoon of Mis on graphene. Light 
(dark) circles indicate C atoms in sublattice A (B). Lines 
show the NN connections, where the electrons hop. Arrows 
illustrate Mis located on top of C atoms, which are either (a) 
distributed equally over both sublattices and AF ordered, or 
(b) located only in sublattice A and FM ordered. 



Magnetic impurities in graphene can either arise 
as effective magnetic moments of non-magnetic 



impurities 



,10.11 



or from intrinsically Mis. The feasi- 



bility of introducing Mis by embedding transition metal 
ions has been studied by ab-initio methodsi^ii^ and can 
be achieved using a scanning tunneling microscope»^ 
Recently, Kondo physics indicating strong coupling be- 
tween conduction electrons and Mis has been reported,— 
suggesting that defects indeed offer a route to tune the 
electronic properties of graphene. Single Mis lead to 
spin-polarized mid-gap statesi^ and have been suggested 
to act as a spin valveJ^ The magnetic (RKKY) coupling 
between impurities is expected to be ferromagnetic (FM) 
[anii-ferromagnetic (AF)] for impurities in the same 
[different] sublattices jiSiii^ii^r^ We are going to discuss 
the ordering of Mis mediated by the electrons and the 
spectral properties of these electrons, which then move 
in a background of ordered Mis. 

The simplest description of graphene is a tight-binding 
Hamiltonian given by the hopping of the carbon- tt^ elec- 
trons between nearest neighbor (NN) sites of a honey- 
comb lattice (see Fig. [T|). The symmetry between the 
two C atoms in the unit cell leads to the Dirac cones 
making up the Fermi surface (FS), located at the points 
A'l and K2 [see inset in Fig. 2(a)| . If this symmetry is 
broken, e.g., by an energy difference between the sublat- 
tices, a gap can be opened, and dispersive edge states 
within the gap were proposed to arise at domain bound- 
aries^i or edges^^ of areas with broken sublattice sym- 
metry. The sublattices could be made inequivalent by a 
substrate with a matching lattice constant^S, or via non- 
magnetic impurities, but only if they mostly occupy sites 
in only one sublattice. ^^ As it will be shown in this letter, 
the situation is different - and possibly simpler - for Mis. 
Due to their AF coupling, they turn out to be oriented in 
opposite directions in the two sublattices and an electron 
that gains magnetic energy in one sublattice, looses it in 
the other. Thus, the sublattices are inequivalent for ei- 
ther spin projection, and both are expected to develop a 
gap regardless of the impurity distribution. If, however, 
Mis or a magnetic substrate can be engineered to inter- 
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FIG. 2: (Color online) Dispersion for electrons that couple 
with J = 1.5 eV to perfectly ordered MI spins, (a) Eq. ([3]) 
for AF ordered Mis on all sites. The curves for up and down 
spins lie on top of each other, (b) Eq. (J4| for Mis on one 
sublattice only; all impurity spins (IS) are parallel and 'up'. 
The dotted line in each panel gives the bands for J = 0. The 
inset in (a) shows the path through the first Brillouin zone. 




act predominantly with only one sublattice, electrons at 
the top (bottom) of the valence (conduction) band will 
be shown to be spin polarized. 

The itinerant electrons interacting with randomly lo- 
cated Mis - represented by an onsite magnetic moment 
Sj^ - are described by the Hamiltonian 



^ = -* E (4-^J.- + ^■'^•) - 2 J E ^i • Si 



(1) 



(i-j),' 



where c[^ creates an electron with spin (T,(i,j) are NN 
sites in the honeycomb lattice, J gives the strength 
of the magnetic interaction between the impurities and 
the electrons, and the spin of the electrons is given by 
Si = X]q /3 Cj Q.'''a,0"^i,^- '^^^ system is half-filled, i.e., 
one electron per site, and t = 3 eV.— i^l To simplify 
the calculations and to avoid the 'sign problem' in the 
Monte Carlo (MC) simulations, we treat the Ml spins as 
classical vectors. This approximation is valid for large 
impurity spins S* > 1, as is expected to hold for transi- 
tion metal ions adsorbed on graphene,— and the local- 
ized spins then simply act as a spin-dependent impurity 
potential on the itinerant electrons. The magnitude of 
the spin is absorbed into the coupling constant J, giving 
2 J sa 3 eV.— Coulomb repulsion between the conduction 
electrons is neglected, as has been found appropriate in 
many physical situations J^i^ Since our lattices are con- 
siderably larger than the distances between impurities, 
finite-size effects are not expected to be severe.— 



II. RESULTS AND DISCUSSION 

For J = 0, the Hamiltonian can be easily diagonalized 
and the band structure along the path K\-T-M-K2 is 
shown as a dotted line in Fig. [2] The two C atoms of the 
basis lead to two states per momentum k, with energies 



±e(k)=±Wl 



-4cos^cos(V3fc,/2)+4cos2:^. (2) 



FIG. 3: (Color online) Spectral density for down electrons, (a) 
IS are in both sublattices; the spectral density for up electrons 
looks identical, (b) IS are in sublattice A only and oriented 
along-f-z, i.e., the down electrons are anti-parallel. j4(k,tj) 
for up electrons is the mirror image about the chemical po- 
tential pi, see Fig. 2(b) The white lines indicate the bands of 
translationally invariant systems described by (a) Eq. ((3]) and 
a renormalized J = cimpJ and (b) Eq. Q and J = 2cimpJ. 



Parameters are: N^ 
J = 1.5 eV = t/2. 



X iVy = 72 X 72, 864 = 72^/6 impurities. 



The momenta k = [k^, ky) allowed for iV^ 



X Ny are k^^ 



Atttix / {Nxa\/o) with Ux — 0,±1, . . . , Nx/2, and ky = 
ATTUy/Ny with Uy = 0, ±1, . . . , Ny/3r^ and the two bands 
are degenerate at the points Ki and K2<^ 

We can also diagonalize Eq. ^ for one MI per lattice 
site with all the Mis in sublattice A (B) parallel to z 
(— z). The spin-degenerate bands are given by 



E{k) = ±VJ2 + e(k)2 



(3) 



Since non-magnetic impurities need to be concentrated 
on one sublattic e to induce a gap, the resulting gap of 
size 2 J [see Fig. 2(a) suggests that Mis might provide a 
simpler way to turn graphene into a semiconductor. 

A different situation arises if all Mi's are in one sublat- 
tice and parallel to z. The energy then depends on the 
spin and for the spin parallel (antiparallel) to z, denoted 
by 'up' ('down'), we find 



J 



EHk)^~-± 



J2 



+ e{ky 



EH-k) = ^± 



J2 

4 



6(k)2 

(4) 



The resulting bands are shown in Fig. 2(b) and not 



surprisingly, the electrons parallel to the localized spins 
are energetically favored. But as one band for each spin 
direction lies below the chemical potential, there is no 
net polarization of the itinerant electrons despite the 
fact that all localized spins are perfectly FM. This fol- 
lows from the degeneracy induced by the two-atom basis 
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FIG. 4: (Color online) Spin structure factor for momenta 
along the path shown in the inset of Fig. 2(a) Symmet- 



ric {S ) and antisymmetric {S~) components indicating FM 
and AF order are shown for the localized IS as well as for the 
itinerant electrons for (a) each sublattice containing half of 
the randomly distributed Mi's and (b) all the Mi's randomly 
distributed in the A sublattice. J — 6 eV, pt — 100, corre- 
sponding to T « 300 K, N^ = Ny ^ 12, Mmp = 24 = 12^/6. 
The data were averaged over (a) six and (b) ten impurity 
configurations. 



and it means that the electrons in the "clean" sublattice 
are polarized in the opposite direction so that the con- 
tributions from the electrons in the two sublattices can- 
cel exactly. While there is thus no ferrimagnetism and 
no net spin polarization of the itinerant electrons, the 
band structure in Fig. 2(b) nevertheless shows that there 



are different energy regions with well developed major- 
ity spins: in particular, the electrons directly below and 
above the Fermi level are polarized. 

As long as J is not too large (J < 2i), the results re- 
main similar even if not all sites of the (sub-) lattice have a 
MI, and the ordered AF leads to dispersive bands instead 
of localized impurity states. To show this, the Hamilto- 
nian given by Eq. ([T]) was diagonalized numerically for 
lattices with N^ x Ny == 72 x 72 = 5184 sites containing 
iVij^p = 864 = 72^/6 randomly distributed Mi's. The 
IS were set into a magnetically ordered state, where all 
spins in sublattice A (B) are parallel (antiparallel) to 
z. The spectral function A{'k,uj) is shown in Fig. 3(a) 
for equal numbers of Mi's in both sublattices. (Notice 
that the bands for spin up and down are here degener- 
ate.) Similar to the case of non-magnetic impuritieSj^ 
the most prominent change is that the gap is renormal- 
ized from 2J to 2ciinpJ for an impurity concentration 
Cimp — Ni^p/{NxNy). This can be seen by comparing 
A(k, w) and the white line given by Eq. ^ with a renor- 
malized J = CjmpJ. The inset shows the resulting gap 
of 2 J = 2cinip^ = 0.5 eV. Doping with Mis thus turns 
graphene into a semiconductor via the AF ordering of 
both the Mis and the graphene electrons. 

In Fig. |3(b)[ we present the band dispersion for the 
'down' electrons (i.e., antiparallel to the IS) for 864 
Mis randomly distributed in one of the sublattices. 
Since the sublattice containing the Mis effectively has 
a doubled impurity concentration of CcS — Scimp — 
Nirapl{NxNy/2), the renormalized gap resulting from 

Eq. dH) is given by J = 2cimpJ' = 0.5 eV, the same as 



for impurities in both sublattices. The gap is thus in- 
dependent of the percentage of impurities found in each 
sublattice. The inset of Fig. |3(b)| reveals that the gap 
now opens above the chemical potential /i, while the gap 
for spin 'up' (not shown) is below //. (The up and down 
bands are mirror images about fj,.) When all Mi's are 
in one sublattice, the highest occupied states are thus 
clearly spin polarized, leading to the spin-polarized half- 
semiconductor described in Ref. 0. 

If the MI distribution is not perfect and a fraction 
p < 1 of the Mis are in one sublattice and 1 — p > 
in the other, the exact degeneracy of up and down states 
at Ki and K2 is hfted and a gap of order 4(1 — p)Jcimp 
opens between the highest occupied 'down' states and the 
lowest empty 'up' states. The gaps for each spin direc- 
tion bands remain asymmetric about /x as long asp 3> 0.5 
(or p <C 0.5), thus preserving the spin polarization for a 
reduced energy window 2(2p— 1) Jcimp. Again, the trans- 
lationally invariant and the dilute systems behave simi- 
larly: The ordered state of the Mis preserves coherent 
electron motion in highly dispersive bands (albeit with a 
finite mass), and we only find localized impurity states 
with small dispersion for very large J > 2t. 

After discussing electrons interacting with perfectly or- 
dered spins, we present results of unbiased MC simula- 
tions. They were performed on the angles 9\ and (pi defin- 
ing the orientation of the IS in the Hamiltonian Eq. ([1]) ; 
the probability of a spin configuration is obtained by di- 
agonalizing the effective electronic Hamiltoniaui^ Stud- 
ies for lattice sizes up to 18 x 18, for 150K< T < 1500K, 
for various values of J, and for different impurity concen- 
trations showed indications of long-range magnetic order. 

Each momentum k has a spin-structure factor SaO<-) 
for sublattice A and SbO^) for sublattice B. They in- 
dicate the magnetic order within the sublattice, with a 
signal at the F-point (k ~ 0) signifying FM order be- 
tween all spins in the same sublattice. Sa and Sb can 
be combined into symmetric (5"*" = Sa + Sb) and an- 
tisymmetric {S^ — Sa — Sb) combinations, which give 
information about the order between the sublattices: If 
the spins in sublattice A and B are parallel (antiparal- 
lel), this is revealed by a signal in 5+ (5"). In Fig. 21 
S~^ and S~ are shown as a function of the momentum. 
One sees clear peaks at F, indicating that the itinerant 
electrons mediate order between the localized spins and 
that all spins within a sublattice are FM. For Mis located 
in both sublattices, the signal is in the S~ channel [see 
Fig. 4(a)| , which shows AF order between the sublattices. 

If sublattice B is free from Mis, it does not con- 
tribute to the impurity-spin structure factor, leading to 
S^ = S^ ^ Sa, and the signal at F indi cates FM order 
within the Ml-doped sublattice [see Fig. |4(b)| . The MC 
simulations thus clearly demonstrate FM order within 
sublattices and AF order between them, justifying the 
choice of ordered ground-state configurations discussed 
abovei^ The structure factor of the itinerant electrons 
develops a peak at F in the AF S~ channel regardless 
of the MI distribution. Since there is no signal in the 



symmetric channel S'^ (it has been verified that this also 
holds for momenta not on the path shown in Fig. U) , fer- 
rimagnetism can be ruled out and the itinerant electrons 
do not develop a net magnetic moment, as discussed pre- 
viously for the perfectly ordered IS. 

For small lattices, the few available momenta mean a 
large finite-size energy spacing of the kinetic energy, and 
an asymmetric gap, as seen in Figs. 2(b) and |3(b)] does 
then not show up very clearly. The smaller lattices acces- 
sible to MC require thus the use of larger values oi J >t 
to allow us to analyze A(k, to). Where the analysis is pos- 
sible, the MC data agree with the ground-state results. 
For Mis in only one sublattice, the spin polarization is re- 
duced by finite-temperature fluctuations, but the states 
around the FS remain predominantly 'down' (w < /i) and 
'up' (a; > /i). The unbiased MC calculations thus corrob- 
orate the results for magnetically ordered configurations. 



III. CONCLUSIONS 

We have discussed itinerant electrons in graphene in- 
teracting with Mis localized on top of the C atoms. We 



found that the electrons always develop AF correlations 
and that a gap opens for both spin projections, regardless 
of the distribution of the Mi's between the two sublat- 
tices. Domain boundaries or edges of nanoribbonsi2ii22 
are expected to support edge states within the gap. For 
Mi's randomly distributed between the sublattices, we 
find a standard semiconductor with a gap that can be 
tuned by adjusting the impurity concentration. While 
the carriers acquire a mass, their dispersion remains al- 
most that of clean graphene. If all or most Mi's are in 
one sublattice^^ the highest occupied (lowest unoccupied) 
states are spin polarized antiparallel (parallel) to the Mis, 
giving rise to a spin-polarized semiconductor that could 
be useful to inject polarized electrons into devices. We 
find that Mis with large spin located on top of C atoms 
are the most likely candidates to produce these effects. 
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